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Abstract 

Background: Gliomas account for almost 80% of primary malignant brain tumors. Epidermal growth factor (EGF) is an 
interesting research candidate in which to look for genetic polymorphisms because of its role in mitogenesis and 
proliferation. Extensive studies have found that a single nucleotide polymorphism (SNP) +61G/A (rs4444903) in the EGF 
gene is associated with the susceptibility of glioma, however, the results have been controversial. Furthermore, the 
association between EGF +61G/A polymorphism with the development and grade progress of glioma has not been 
established. 

Methods: We examined the association of EGF +61G/A polymorphism and glioma by performing a meta-analysis. Nine 
studies testing the associations between EGF +61G/A polymorphism and risk of glioma with 1758 cases and 2823 controls 
were retrieved. Odds ratios (ORs) and 95% confidence intervals (CIs) were used to assess the strength of the association. The 
pooled ORs were performed for the allele model, codominant model, dominant model, and recessive model, respectively. 

Results: Overall, this meta-analysis showed significant associations between the EGF +61G/A polymorphism and glioma 
susceptibility in all four genetic models. However, in the stratified analysis by the grade of glioma, we only found this 
association existed in patients with Grade IV glioblastoma, but not in patients with Grade l-lll glioma. We further compared 
EGF +61G/A polymorphism in patients with glioblastoma and Grade l-lll glioma accordingly, the stronger association 
between the EGF +61G/A polymorphism and the malignancy of glioma was found. 

Conclusions: The results of this meta-analysis suggested that the EGF +61G/A polymorphism is associated with both the 
susceptibility of glioma and the malignance of glioma. 
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Introduction 

Gliomas are the most common primary tumors in the central 
nervous system [1]. Gliomas have complex cellular composition 
and inherent resistance to the treatments. Carcinogenesis of 
glioma is a complex pathological process, in which environmental 
and genetic factors interact with each other [2] . Ionizing radiation 
is a known factor which is clearly associated with brain tumors 
development risk [3]. However, many patients without known 
environmental factors eventually developed glioma. Therefore, 
genetic predisposition may contribute to the process of glioma 
carcinogenesis. According to the World Health Organization 
(WHO) criteria, gliomas are histologically classified into Grades I 
to IV [4] . The higher the glioma grade is, the more resistant to the 
treatments and the shorter median survival time the patient will 
have [5] . The genetic heterogeneity of glioma also accounts for the 
current limitations in predicting patient prognosis based on 



histologic classification of glioma grade alone [6,7] and suggests 
that classification of gliomas from the same grade according to 
their genetic phenotype will provide a more accurate prediction of 
both patient survival time and the response to the therapy [8,9] . 

The epidermal growth factor (EGF) gene, which is a member of 
the EGF superfamily, has been demonstrated to activate cell 
proliferation and stimulate mitogenesis in epidermal tissues [10]. It 
is a natural focus of interest in cancer studies because of its ability 
to influence mitogenesis and proliferation, and its overexpression is 
commonly seen in human cancers, such as glioma, pancreatic 
cancer, breast cancer and hepatocarcinoma [11,12]. Some studies 
reported that EGF plays a critical role in gliomagenesis by 
influencing cell proliferation, inhibiting apoptosis, and inducing 
cell differentiation, via binding to its specific receptor, the EGF 
receptor (EGFR) [13,14]. Previously, Shahbazi et al. identified a 
functional single nucleotide polymorphism involving guanine (G) 
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Figure 1. Flow chart of the study selection process. 

doi:10.1371/journal.pone.0095139.g001 

to adenine (A) substitution at position 61 in the 5 '-untranslated 
region of the EGF gene (rs4444903), and the substitution increased 
production of EGF in cultured peripheral blood mononuclear cells 
[15]. Thus, this genetic polymorphism may contribute to 
interindividual differences of EGF expression and subsequently 
the tumor predisposition and aggressiveness. 



To date, several studies have found the association between the 
EGF +61G/A polymorphism with the susceptibility of glioma, 
however, these studies showed inconsistent results [16-19]. Some 
molecular epidemiological studies suggested that the inconsistency 
might be caused by the diverse ethnicity of patients [16,17], 
however, we hypothesize that the diverse distribution of EGF gene 
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frequencies among the different ethnic groups may cause the 
inconsistency [20,21]. In addition, previous researches did not 
analyze the association in different glioma grade. The compre- 
hensive analysis of glioma grade and genetic phenotype may more 
accurately predict the prognosis of patients. Therefore, a meta- 
analysis of nine case-control studies involving 1,758 cases and 
2,823 controls was performed to drive a more precise estimation of 
the association of EGF +61G/A polymorphism with susceptibility 
to and severity of the glioma. 

Materials and Methods 

The procedures of this meta-analysis adhered to the PRISMA 
Statement guidelines [22] (Supporting Information Checklist SI). 

Literature search strategy for identification of the studies 

We did a literature search in online database such as PubMed, 
Embase and the Chinese knowledge library using the search terms 
"epidermal growth factor", "polymorphism", and "glioma", along 
with additional terms such as "EGF, SNP, Glioblastoma, 
astrocytoma, oligodendroglioma", and all possible combinations 
(from January 1, 2000 to September 5, 2013). The search was 
limited to English-language articles. In order to search more 
deeply, we reviewed the references of the selected articles to 
retrieve data that we could have ignored in the initial search. Of 
the studies with overlapping data published by the same 
researchers, we selected the most recent ones with the largest 
number of subjects. 

Inclusion and exclusion criteria 

The studies for the meta-analysis were selected when they 
satisfied the following criteria: studies about EGF +61G/A 
polymorphism and glioma risk, case-control studies in humans, 
studies containing useful genotype frequencies or odds ratio (OR) 
data, studies containing classification of glioma grade. Major 
exclusion criteria were (a) no control population, (b) no available 
genotype frequency, (c) no explicit glioma grade classification, and 
(d) duplication of the previous publication. 

Data extraction 

Information was extracted from all eligible articles indepen- 
dently by two of the authors following the criteria listed above (Xin 
Chen and Guang Yang). Disagreement was discussed and resolved 
between the two authors. 

For each study, the following data were considered: the first 
author's last name, year of publication, country of origin, ethnicity 
and numbers of genotyped cases and controls. We did not define 
any minimum number of patients to include a study in our meta- 
analysis. In the event that a study presented subpopulations, these 
were taken to be different studies. 

Statistical analysis 

The estimate of association between the EGF +61G/A 
polymorphism and glioma risk was assessed using the fixed-effect 
method [23] which calculates the crude ORs and the correspond- 
ing 95% CIs for individual studies and the global association. In 
this meta-analysis, small letter "a" was defined as the minor allele 
or the mutant allele of EGF +61G/A polymorphism, while capital 
letter "A" was defined as the major allele or the common allele of 
EGF +61G/A polymorphism. The pooled ORs were performed 
for the allele model (a versus A), the codominant model (aa versus 
AA; Aa versus AA), the dominant model (aa/Aa versus AA), and 
the recessive model (aa versus Aa/AA), respectively. Stratified 
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analyses were performed by tumor grade. When the test was 
heterogeneous, the random-effect method [24] was applied. 

The analysis of heterogeneity between studies was performed by 
the Q statistic, with p-values<0.05 indicating significant hetero- 
geneity [25]. Sensitivity was performed by omitting individual 
studies in order to assess the stability of results. Begg's funnel plots 
were used to detect publication bias. Additionally, Egger's linear 
regression test which measures funnel plot asymmetry was also 
used to evaluate the publication bias, If P<0.05, statistically 
significant publication bias existed [26]. All analyses were 
calculated using the STATA Version 11.0 software (version 11; 
Stata Corporation, College Station, Texas). 

Results 

Characteristics of included Studies 

Figure 1 outlines our study selection process. Briefly, a total of 
88 potentially relevant articles were systematically identified after 
an initial search. After reading the titles and abstracts, we removed 
58 articles not examining the association between EGF polymor- 
phism and glioma risk, 11 articles not examining EGF +61G/A 
polymorphism, and 9 articles belong to review and meta-analysis. 
After reading full texts of the remaining 1 0 articles regarding the 
association between the rs4444903 polymorphism on the EGF 
gene and glioma, one of these was excluded for not providing allele 
frequencies in glioma grade subgroup needed for OR calculation 
[27]. Finally, a total of 9 case— control studies from 5 countries 
involving 1,758 cases and 2,823 controls were accepted for the 
pooled analysis. Table 1 presents the main characteristics of these 
studies. 

Among the 9 publications, there were five studies of Caucasian 
descendants [20,28-31] and four studies of Asian descendants 
[21,32-34]. Gliomas were confirmed histologically or pathologi- 
cally in all included studies. A classic polymerase chain reaction- 
restriction fragment length polymorphism (PCR-RFLP) assay was 
performed in most studies [21,29-32,34], only one study used 
polymerase chain reaction-ligation detection reaction (PCR-LDR) 
method [33], one study used quantitative reverse transcription- 
PCR (qRT-PCR) technique [28], and one study used Real-time 
PCR method [20]. The controls were mainly hospital populations 
and all studies indicated that the distribution of genotypes in the 
controls were consistent with Hardy- Weinberg equilibrium (HWE) 
except for two studies [20,34] . 

Association between EGF +61G/A polymorphism and 
glioma risk 

The evaluation of association between EGF +61G/A polymor- 
phism and glioma risk are presented in Table 2. There were no 
significant heterogeneity between studies in all comparisons. 
Overall, significant associations between the EGF +61G/A 
polymorphism and glioma risk were observed in all genetic 
models (a versus A, OR = 1.22, 95% CI = 1 . 1 1-1.33; aa/Aa versus 
AA, OR=1.21, 95% CI= 1.06-1.38) (Fig. 2 and Fig. SI). 
However, in the stratified analysis by the grade of glioma, we 
found the EGF +61G/A polymorphism was significantly associ- 
ated with increased risk of glioma among patients with glioblas- 
toma (Grade IV) (a versus A, OR= 1.31, 95% CI= 1.11-1.55; aa/ 
Aa versus AA, OR=1.31, 95% CI = 1.1 1-1.55), and but not 
among patients with Grade I-III glioma (Grade I: a versus A, 
OR= 1.20, 95% CI = 0.91-1.59; aa/Aa versus AA, OR= 1.18, 
95% CI = 0.81-1.74; Grade II: a versus A, OR= 1.15, 95% 
CI = 0.98-1.35; aa/Aa versus AA, OR= 1.16, 95% CI = 0.94- 
1.43; Grade III: a versus A, OR= 1.08, 95% CI = 0.89-1.31; aa/ 
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Figure 2. Forest plot of glioma risk for EGF +61 G/A polymorphism (a versus A). The squares and horizontal lines correspond to the study- 
specific OR and 95% CI. The area of the squares reflects the study-specific weight (inverse of the variance). The diamond represents the pooled OR 
and 95% CI. 

doi:1 0.1 371 /journal.pone.00951 39.g002 



Aa versus AA, OR= 1.14, 95% CI = 0.88-1.48) (Fig. 3 and Fig. 
S2). 

At the same time, on the basis of glioma grade subgroup, we 
further compared the EGF +61 G/A polymorphism in patients 
with glioblastoma and Grade I-III glioma, and found the OR 
value was greater than 1 (a versus A, OR= 1.28, 95% CI= 1.08— 
1.51; aa/Aa versus AA, OR= 1.30, 95% CI= 1.03-1.63), in the 
sense that by contrast to Grade I-III glioma, glioblastoma has the 
stronger associations with EGF +6 1 G/ A polymorphism (Fig. 4 and 
Fig. S3). 

Sensitivity analysis 

To assess the influence of each individual study on the pooled 
ORs, we performed a sensitivity analysis by omitting individual 
study. The analysis results suggested that no individual study 
significantly affected the pooled ORs under the allele model (a 
versus A) of EGF +61 G/A polymorphism (Fig. 5), which indicates 
that our results are statistically robust. 

Publication bias 

Begg's funnel plot and Egger's test were performed to assess the 
publication bias of included literatures. The shapes of the funnel 



plots in all genetic models of EGF +6 1 G/ A polymorphism seemed 
symmetrical (Fig. 6), suggesting the lack of publication bias in the 
meta-analysis. Then, the Egger's test was used to provide statistical 
evidence of funnel plot symmetry. As expected, the results shown a 
significantly statistical evidence of no publication bias (t = —0.47, 
P = 0.656 for a versus A). 

Discussion 

Cell growth dysregulation is one important characteristic of 
tumors. EGF/EGFR signaling dysregulation often occurs in 
gliomas and plays a key role in gliomagenesis [35,36]. EGF can 
activate a cascade of intracellular signaling molecules, these 
molecules are important for cell proliferation, survival, migration, 
and differentiation [37]. The abnormalities of EGF/EGFR- 
mediated pathway are often associated with poor clinical outcome 
[38]. Therefore, it is logical to hypothesize that functional genetic 
polymorphisms involving EGF or EGFR may predispose to glioma 
development. 

Many retrospective studies have shown that polymorphism in 
EGF or EGFR genes strongly correlates with susceptibility to 
glioma. Some studies have reported the association between 
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Figure 3. Forest plot of glioma risk for EGF +61G/A polymorphism (a versus A) by glioma grade. The dots and horizontal lines 
correspond to the study-specific OR and 95% CI. The diamond represents the pooled OR and 95% CI. 
doi:1 0.1 371 /journal.pone.00951 39.g003 



polymorphism in EGF +61G/A and glioma susceptibility 
[20,21,28-34], but the results of these studies are not consistent. 
Various teams considered the inconformity is derived from 
different ethnicity [16,17]. However, we hypothesize that this 
racial discrepancy could attribute to the diverse distribution of 
EGF gene frequencies among the different ethnic groups. In this 
meta-analysis, on the basis of collecting more studies than previous 
systemic review, we tried to eliminate the racial difference by 
defining the small letter "a" as the minor allele or the mutant allele 
of EGF +61G/A polymorphism, the capital letter "A" as the 
major allele or the common allele of EGF +61G/A polymorphism. 
As we expected, the analysing results showed the consistent and 
significant associations between the EGF +61G/A polymorphism 
and glioma susceptibility in in all genetic models. 



From the above, we concluded that EGF +61G/A polymor- 
phism is associated with the susceptibility to glioma. As the 
literature reported, EGF +61G/A polymorphism is also probably 
associated with the survival time of glioma patients [30]. Since 
glioma grade is the most informative factor for stratification into 
subgroups with different prognoses [5], We speculated that EGF + 
61G/A polymorphism may be associated with the glioma grade. 
Becuase many studies did not estimate the association between 
EGF +61G/A polymorphism and the risk of glioma at different 
grade. So a strict analysis of the association of EGF +61G/A 
polymorphism and glioma risk should be performed. 

Our approach allowed us to look for potential glioma grade 
differences in the association between EGF +61G/A polymor- 
phism and glioma risk. Analysis of glioma grade subgroups showed 
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Figure 4. Forest plot of EGF +61G/A polymorphism associated with glioblastoma and Grade I — ■ ■ ■ glioma risk (a versus A). The squares 
and horizontal lines correspond to the study-specific OR and 95% CI. The area of the squares reflects the study-specific weight (inverse of the 
variance). The diamond represents the pooled OR and 95% CI. 
doi:10.1371/journal.pone.0095139.g004 



Meta-analysis fixed-effects estimates (exponential form) 
Study ommited 



Bao(2010) 
Wang(2010) 
Liu(2009) 
Li(2012) 
Pinto(2009) 
Costa(2007) 
Silveira(2012) 
Bhowmick(2004) 
Vauleon(2007) 



1.05 



1.11 



1.22 



1.33 



1.38 



Figure 5. Sensitivity analysis of the association between EGF +61G/A polymorphism and glioma risk (a versus A). Results were 
computed by omitting each study in turn. Meta-analysis fixed-effects estimates (exponential form) were used. Open circle indicates the pooled OR, 
given named study is omitted. The two ends of the dotted lines represent the 95% CI. 
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a strong association between subjects carrying at least one "a" 
allele with the risk of glioblastoma (Grade IV), but this association 
was not found in Grade I-III glioma. We further compared the 
EGF +61 G/A polymorphism in patients with glioblastoma and 
Grade I-III glioma, and found that glioblastoma has a stronger 
association with EGF +61 G/A polymorphism compared with 
Grade I-III glioma, this means there is a significant association 
between the EGF +61 G/A polymorphism and the malignancy of 
glioma. 

We think this result will strengthen the previous knowledge 
about EGF +61 G/A polymorphism and glioma risk. This study 
has an important implication for glioma progress prediction that 
genotyping for the EGF +61 G/A variant might help in identifying 
patients at high risk of glioblastoma. However, it needs to be 
emphasized that our meta-analysis involved only one study with 
Caucasian Grade I-III glioma population, so this association should 
be re-evaluated in studies with larger sample sizes. Therefore, the 
results in the present study should be interpreted with caution. 

There are some limitations that should be considered to 
interpret this result. Confounding factors may affect the meta- 
analysis. Firstly, detailed information such as age and gender in 
cases and controls of different genotypes in some studies were not 
available, which limited further analysis. Secondly, the numbers of 
published studies were not enough for a comprehensive analysis, 
particularly for different grades of glioma, we could not obtain 
enough statistical power to investigate the real association. Because 
of the limited analysing data, we can only conclude that EGF + 
61 G/A polymorphism is related to high grade gliomas. In current 
results, the correlativity between EGF +61 G/A polymorphism and 
low grade glioma is not clear. Until now, little is known about its 
function in primary and secondary glioblastomas. As Kunihiko 
Watanabe reported that immunoreactivity for the EGFR prevailed 
in primary glioblastomas but was rare in secondary glioblastomas 
[39], it infers that disturbance of EGF-EGFR signalling pathway 
including SNP in EGF and EGFR genes may generally participate 
in primary glioblastomas instead of the secondary glioblastomas. 
This speculation corresponds with our conclusion more or less, 
however, whether EGF +61 G/A polymorphism mutates during 
the development from low grade to high grade gliomas needs 



further experimental explorations, more studies are needed to 
detect the EGF +61 G/A polymorphism and their associations with 
those classified gliomas. 

Despite of these limitations, our meta-analysis still has clear 
advantages over previous studies. To the best of our knowledge, 
this is the first meta-analysis of the association between EGF + 
61 G/A gene polymorphisms and the susceptibility to glioma 
subject to different grades. In addition, data were abstracted from 
different studies, which significandy increased statistic power of 
this meta-analysis. Meanwhile, studies included in this meta- 
analysis were case-control studies and contained available 
genotype frequency, which met our selection criterions. 

In summary, this meta-analysis suggested that the EGF +61G/ 
A polymorphism is associated with glioma development risk. 
Additionally, we found that this association was more prominent 
in glioblastoma than in Grade I-III glioma, showing that the EGF 
+61 G/A polymorphism is not only associated with the suscepti- 
bility to glioma, but also with the malignance of glioma. Future 
well designed larger studies are warranted to validate this 
association in different grade populations of glioma. 

Supporting Information 

Figure SI Forest plot of glioma risk for EGF +61G/A 
polymorphism (aa/Aa versus AA). The squares and 
horizontal lines correspond to the study-specific OR and 95% 
CI. The area of the squares reflects the study-specific weight 
(inverse of the variance). The diamond represents the pooled OR 
and 95% CI. 
(TIF) 

Figure S2 Forest plot of glioma risk for EGF +61G/A 
polymorphism (aa/Aa versus AA) by glioma grade. The 

dots and horizontal lines correspond to the study-specific OR and 
95% CI. The diamond represents the pooled OR and 95% CI. 
(TIF) 

Figure S3 Forest plot of EGF +61 G/A polymorphism 
associated with glioblastoma and Grade I-III glioma 
risk (aa/Aa versus AA). The squares and horizontal lines 
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correspond to the study-specific OR and 95% CI. The area of the 
squares reflects the study-specific weight (inverse of the variance). 
The diamond represents the pooled OR and 95% CI. 
(TIF) 

Checklist SI PRISMA Checklist. 

(DOC) 
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